Unitarity Bounds for Gauged Axionic Interactions and the 

Green-Schwarz Mechanism 



'Claudio Coriano "Marco Guzzi and a Simone Morelli 



a Dipartimento di Fisica, Universita del Salento 
and INFN Sezione di Lecce, Via Arnesano 73100 Lecce, Italy 

b Department of Physics and Institute of Plasma Physics 
University of Crete, 71003 Heraklion, Greece 

Abstract 

We analyze the effective actions of anomalous models in which a four-dimensional ver- 
sion of the Green-Schwarz mechanism is invoked for the cancellation of the anomalies, and 
we compare it with those models in which gauge invariance is restored by the presence 
of a Wess-Zumino term. Some issues concerning an apparent violation of unitarity of 
the mechanism, which requires Dolgov-Zakharov poles, are carefully examined, using a 
class of amplitudes studied in the past by Bouchiat-Iliopoulos- Meyer (BIM), and elabo- 
rating on previous studies. In the Wess-Zumino case we determine explicitly the unitarity 
bound using a realistic model of intersecting branes (the Madrid model) by studying the 
corresponding BIM amplitudes. This is shown to depend significantly on the Stiickelberg 
mass and on the coupling of the extra anomalous gauge bosons and allows one to identify 
Standard-Model-like regions (which are anomaly-free) from regions where the growth of 
certain amplitudes is dominated by the anomaly, separated by an inflection point which 
could be studied at the LHC. The bound can even be around 5-10 TeV's for a Z' mass 
around 1 TeV and varies sensitively with the anomalous coupling. The results for the WZ 
case are quite general and apply to all the models in which an axion-like interaction is 
introduced as a generalization of the Peccei-Quinn mechanism, with a gauged axion . 
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1 Introduction 



The cancellation of gauge anomalies in the Standard Model (SM) is a landmark of modern 
particle theory that has contributed to shape our knowledge on the fermion spectrum, its chiral 
charges and couplings. Other mechanisms of cancellation, based on the introduction of both 
local and non-local counterterms, have also received a lot of attention in the last two decades, 
from the introduction of the Wess-Zumino term in gauge theories [T] (which is local) to the 
Green-Schwarz mechanism of string theory [2] (which is non-local). The field theory realization 
of this second mechanism is rather puzzling also on phenomenological grounds since it requires, 
in four dimensions, the non-local exchange of a pseudoscalar to restore gauge invariance in the 
anomalous vertices. In higher dimensions, for instance in 10 dimensions, the violation of the 
Ward identities due to the hexagon diagram is canceled by the exchange of a 2- form [21 [3] . In 
this work we are going to analyze the similarities between the two approaches and emphasize 
the differences as well. We will try, along the way, to point out those unclear aspects of the 
field theory realization of this mechanism - in the absence of supersymmetry and gravitational 
interactions - which, apparently, suffers from the presence of an analytic structure in the energy 
plane that is in apparent disagreement with unitarity. Moving to the WZ case, here we show 
that the restoration of gauge invariance in the corresponding one-loop effective Lagrangian via 
a local axion counterterm is not able to guarantee unitarity beyond a certain scale, although 
this deficiency is expected [U |5], given the local nature of the counterterm. In the GS case, 
the restoration of the Ward identities suffers from the presence of unphysical massless poles 
in the trilinear gauge vertices that, as we are going to show, are similar to those present in 
a non-local version of axial electrodynamics, which has been studied extensively in the past 
[6] with negative conclusions concerning its unitarity properties. In particular, in the case of 
scalar potentials that include Higgs- axion mixing, the phenomenological interpretation of the 
GS mechanism remains problematic in the field theoretical construction. 

We comment on the relation between the two mechanisms, when the axion is integrated out 
of the partition function of the anomalous theory, and on other issues of the gauge dependence 
of the perturbative expansions, which emerge in the different formulations. In the second part 
of this work we apply our analysis to a realistic model characterizing numerically the bounds 
in effective actions of WZ type and discuss the possibility to constrain brane and axion-like 
models at the LHC. 
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1.1 WZ and GS counterterms 

Anomalous abelian models are variations of the standard model in which the gauge structure 
of this is enlarged by one or more abelian factors. The corresponding anomalies are canceled by 
the introduction of a pseudoscalar, an axion (b), that couples to 4- forms FjAFj (via b/MFjAFj, 
the Wess-Zumino term) of the gauge fields (/, J) that appear both in ordinary (I = J) and 
mixed (/ 7^ J) anomalies. M is a scale that is apparently unrelated to the rest of the theory and 
simply describes the range in which the anomalous model can be used as a good approximation 
to the underlying complete theory. The latter can be resolved at an energy E > M, by using 
either a renormalizable Lagrangian with an anomaly-free chiral fermion spectrum or a string 
theory. The motivations for introducing such models are several, ranging from the study of 
the flavor sector, where several attempts have been performed in the last decade to reproduce 
the neutrino mixing matrix using theories of this type, to effective string models, in which the 
extra U(l) abound. We also recall that in effective string models and in models characterized 
by extra dimensions the axion (b) appears together with a mixing to the anomalous gauge 
boson (dbB), which is, by coincidence, natural in a (Higgs) theory in a broken phase. In a way, 
theories of this type have several completions at higher energy [5]. 

Coming to the specific models that we analyze, these are complete MLSOM-like [3, [8] models 
with three anomalous U(l) P, \10\, while most of the unitarity issues are easier to address in 
simple models with two U(l) [TTj . In our phenomenological analysis, which concerns only 
effective actions of WZ type, we will choose the charge assignments and the construction of 
[12J, but we will work in the region of parameter space where only the lowest Stiickelberg mass 
eigenvalue is taken into account, while the remaining two extra Z' decouple. This configuration 
is not the most general but is enough to clarify the key physical properties of these models. 

2 Anomaly cancellation and gauge dependences: the GS 
and the WZ mechanisms in field theory 

In this section we start our discussion of the unitarity properties of the GS and WZ mechanisms, 
illustrating the critical issues. We illustrate a pure diagrammatic construction of the WZ 
effective action using a set of basic local counterterms and show how a certain class of amplitudes 
have an anomalous behavior that grows beyond their unitarity limit at high energy. The 
arguments being rather subtle, we have decided to illustrate the construction of the effective 
action for both mechanisms in parallel. A re-arrangement of the same basic counterterms of the 
WZ case generates the GS effective action, which, however, is non-local. The two mechanisms 
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are different even if they share a common origin. In a following section we will integrate out the 
axion of the WZ formulation to generate a non-local form of the same mechanism that resembles 
more closely the GS counterterm. The two differ by a set of extra non-local interactions in their 
respective effective actions. One could go the other way and formulate the GS mechanism in 
a local form using (two or more) extra auxiliary fields. These points are relevant in order to 
understand the connection between the two ways to cancel the anomaly. 

2.1 The Lagrangian 

Specifically, the toy model that we consider has a single fermion with a vector-like interaction 
with the gauge field A and a purely axial- vector interaction with B. The Lagrangian is given 
by 

A) = ~Fl-^ + ^(d,b + M 1 B^ 2 + ^iY(d^ + ieA fi + ig Bl 5 B^, (1) 

where, for simplicity, we have taken all the charges to be unitary, and we have allowed for a 
Stiickelberg term for B, with Mi being the Stiickelberg mass Q- A is massless and takes the 
role of a photon. The Lagrangian has a Stuckelberg-like symmetry with b — > 6 — Mi$b under a 
gauge transformation of B^, 5B^ = <9 m #b. The axion is a singlet under gauge transformations 
of A. We call this simplified theory the "A — B" model. We are allowed not to perform any 
gauge fixing on B and keep the coupling of the longitudinal component of B to the axion, dBb, 
as an interaction vertex. If we remove A, we call the simplified model the "B model". We 
will be interchanging between these two models for illustrative purposes and to underline the 
essential features of theories of this type. 

In the A — B model, the U(1)a gauge freedom can be gauge-fixed in a generic Lorenz gauge, 
with polarization vectors that carry a dependence on the gauge parameter £4, but A being non- 
anomalous we will assume trivially the validity of the Ward identities on vector-like currents. 
This will erase any dependence on £4 both of the polarization vectors of A and of the propagators 
of the same gauge boson. At the same time Chern-Simons (CS) interactions such as AB A Fb 
or AB A Fa, which are present if we define triangle diagrams with a symmetric distribution 
of the partial anomalies of each vertex both in the AVV (axial- vector /vector/vector) and 
AAA cases [TiT [TO] , can be absorbed by a re-distribution of the anomaly. For instance, if we 

1 Even if fl} is not the most general invariant Lagrangian under the gauge group U(1)a x U(1)b, our 
considerations are the same. In fact, since b shifts only under a gauge variation of the anomalous U(X) gauge 
field B (and not under A), the gauge invariance of the effective action under a gauge transformation of the 
gauge field A requires that there are no terms of the type 6-Fa A Fb ■ 
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assume vector Ward identities on the A current and move the whole anomaly to the axial- 
vector currents, then the CS terms can be omitted. The anomalous corrections in the one-loop 
effective action are due to triangle diagrams of the form BAA (AW, with conserved vector 
currents) and BBB (AAA with a symmetric distribution of the anomalies) which require 
two WZ counterterms, given in Swz below, for anomaly cancellation. Since the analysis of 
anomalous gauge theories containing WZ terms has been the subject of various analyses with 
radically different conclusions regarding the issue of unitarity of these theories, we refer to 
the original literature for more details [H [131 [HI [15]. Our goal here is to simply stress the 
relevance of these previous analyses in order to understand the difference between the WZ 
and GS cancellation mechanism and clarify that Higgs- axion mixing does not find a suitable 
description within the standard formulation of the GS mechanism. 

3 Local and non-local formulations 

The GS mechanism is closely related to the WZ mechanism pQ. The latter, in this case, consists 
in restoring gauge invariance of an anomalous theory by introducing a shifting pseudoscalar, an 
axion, that couples to the divergence of an anomalous current. It can be formulated starting 
from a massive abelian theory and performing a field- enlarging transformation pj] so as to 
generate a complete gauge invariant model in which the usual abelian symmetry is accompanied 
by a shifting axion. The original Lagrangian of the massive gauge theory is interpreted as 
the gauge-fixed case of the field-enlarged Lagrangian. The gauge variation of the anomalous 
effective action is compensated by the WZ term, so as to have a gauge invariant formulation 
of the model. We will show next that a theory built in this way has a unitarity bound that we 
will be able to quantify. The appearance of the axion in these theories seems to be an artifact, 
since the presence of a symmetry allows one to set the axion to vanish, choosing a unitary 
gauge. In brane models, in the presence of a suitable scalar potential, the axion ceases to be 
a gauge artifact and cannot be gauged away, as shown in [7]. This point is rather important, 
since it shows that the GS counterterm is unable to describe Higgs- axion mixing, which takes 
place when the Peccei-Quinn [T6l [T7] symmetry of the scalar potential is broken. The reason is 
quite obvious: the GS virtual axion is a massless exchange whose presence is just to guarantee 
the decoupling of the longitudinal component of the gauge boson from the anomaly and that 
does not describe a physical state. But before coming to a careful analysis of this point, let us 
discuss the counterterms of the Lagrangian. 

In the A — B model the WZ counterterms are 
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Figure 2: A typical Bouchiat-Iliopoulos-Meyer amplitude and the axion counterterm to restore gauge 
invariance in the gauge in the WZ effective action. 



which are fixed by the condition of gauge invariance of the Lagrangian. The best way to proceed 
in the analysis of this theory is to work in the gauge in order to remove the B-b mixing |TT] . 
Alternatively, we are entitled to keep the mixing and perform a perturbative expansion of the 
model using the Proca propagator for the massive gauge boson, and treat the bdB term as a 
bilinear vertex. This second approach can be the source of some confusion, since one could be 
misled and identify the perturbative expansion obtained by using the WZ theory with that of 
the GS mechanism, which involves, at a field theory level, only a re-definition of the trilinear 
fermionic vertex with a pole-like counterterm. In the WZ effective action, treated with the 
b — B mixing (and not in the gauge), similar counterterms appear. 

We show in Fig. [1] the vertices of the effective action in the R^ gauge approach, and we 
combine them to describe the process AA — > AA, as shown in FigJ21 Graph a) of Figj2] is a 
typical BIM amplitude [18], first studied in '72 by Bouchiat, Iliopoulos and Meyer to analyze 
the gauge independence of anomaly-mediated processes in the Standard Model. The gauge 
independence of this process is a necessary condition in a gauge theory in order to have a 
consistent S-matrix free of spurious singularities [TT], but is not sufficient to guarantee the 
absence of a unitarity bound. Typically, gauge cancellations help to identify the correct power 
counting (in 1/M\ and in the coupling constants) of the theory and are essential to establish the 
overall correctness of the perturbative computations using the vertices of the effective action. 
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Figure 3: Perturbative expansion of the AA^AA amplitude in the presence of B — b mixing. 

In our example, this can be established as follows: diagram b) of Fig[2] cancels the gauge 
dependence of diagram a) but leaves an overall remnant, which is the contribution of diagram 
a) computed in the unitary gauge (£ — > oo) in which the propagator takes the Proca form 

- -whif (>' - %) ■ (3) 

Due to this cancellation, the total contribution of the two diagrams is 

A x / v (k,p 1 ,p 2 )D xx 'A x ;^'(k,k 1 ,k 2 ), (4) 

where A^ is given by the Dolgov-Zakharov parameterization 

A A (k,p 1 ,p 2 ) = A 6 (s)(k 1 + k 2 ) x e [h, k 2 , u,/j] , (5) 

where the coefficient A 6 (s) in the massless case is A 6 (s) = l/2(7r 2 s). 

We call Ab the triangle with a symmetric distribution of the anomaly (a n /3 for each vertex), 
which is obtained from A^ by the addition of suitable CS terms [HI HO]- The bad behavior of 
this amplitude at high energy is then trivially given by 

with a n = i/(2n 2 ). Squaring the amplitude, the corresponding cross section grows linearly with 
s = k 2 , which signals the breaking of unitarity, as expected in Proca theory, if the corresponding 
Ward identities are violated. A similar result holds for the BBB case. In the alternative 
formulation, in which the b — B term is treated as a vertex, the perturbative expansion is 
formulated diagrammatically as in Fig. [3J Though the expansion is less transparent in this 
case, it is still expected to reproduce the results of the gauge and of the unitary gauge. 
Notice that the expansion seems to generate the specific GS counterterms (Graph 3b)) that 
limits the interaction of the gauge field with the anomaly to its transverse component, together 
with some extra graphs, which are clearly not absorbed by a re-definition of the gauge vertex. 

2 We will use the coincise notation e[X,p, k, v] = t Xa ^ u p a kp and so on. 
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3.1 Integrating out the Stiickelberg in the WZ case 

We can make a forward step and try to integrate out the axion from the partition function and 
obtain the non-local version of the WZ effective action. Notice that this is straightforward only 
in the case in which Higgs- axion mixing is absent. The partition function in this case is given 
by 

Z = J D^D^jDADBDbexp(i(£(tP,iP,A,B,b))) , (7) 
where ( ) denote integration over x and 



with C and Cwz given in ([I]) and (E 
partial integration on the axion gives l 



C — C + Cwz, (8) 
, respectively. Indicating with C b the b sector of C, a 



where 



C b = -~bDb + bJ, (9) 
J = MdB - ^F A A Fa — ^F B A F B , (10) 



and performing the path integration over b we obtain 



j Dbexp (i(C)) = det (□ - M x 2 ) 1/2 exp (~ JD' 1 J^j , (11) 

where 

(JD-^wz = ((.MJB - ^-F A AF A - ^-F B A F B ^j D' 1 (m^B - ^-F A A F A — ^-F B A F B 

1 (12) 

The additional contributions to the effective action are now non-local and are represented by 
the set of diagrams in FigJH Among these diagrams there are two GS counterterms (diagrams 
c) and d)), but there are also other contributions. To generate only the GS counterterms one 
needs an additional pseudoscalar called a in order to enforce the cancellation of the extra terms. 
There are various ways of doing this [T31 El [12] ■ In [H] the non-local counterterm dB\3~ 1 F AF 
of axial QED, which corresponds to the diagrams b) and c) of FiglSl is obtained by performing 
the functional integral over a and b of the following action [19] 



C = ^(i jd + e <B lh )i)-~F 2 B + — ^— F B AF B (a + b) 

4 4Stt z Mi 

+~ (d„b - M x Btf - \ (c> - M x B„f . (13) 



3 We have re-defined the coefficients in front of the counterterms absorbing the multiplicity factors. 
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A A B B A 

1 b rO 1 b ^ B b rO 

S b B 

(cfj B B (e) A (f) 
Figure 4: The structure of the WZ effective action having integrated out the axion. 

a, b ^ a, a, b a, b 

+ WV*- ^ + AAA, + WV" AAA/ 








(a) ffcj fcj fdj 

Figure 5: Effective action with two pseudoscalars, one of them ghostlike, a. 

The integral on a and 6 are gaussians and their contributions to the effective action, after 
integrating them out, are shown in Figj5j Notice that b has a positive kinetic term and a is 
ghost-like. The role of the two pseudoscalars is to cancel the contributions in Fig. [5] a) and d), 
leaving only the contribution given by graphs b) and c), which has the pole structure typical of 
the GS non-local counterterm. Due to these cancellations, the effective action now reduces to 

(dB(x)D- 1 (x-y)F(y)AF(y)) (14) 

besides the anomaly vertex and is represented by interactions of the form b) and c) of Fig. 
El This shows that the WZ and the GS effective actions organize the perturbative expansions 
in a rather different way. It is also quite immediate that the cleanest way to analyze the 
expansion is to use the gauge, as we have already stressed. It is then also quite clear 
that in the WZ case we require the gauge invariance of the Lagrangian but not of the trilinear 
gauge interactions, while in the GS case, which is realized via Eq. ( 1T31) or, analogously, by 
the Lagrangians proposed in [T3l IT4] , it is the trilinear vertex that is rendered gauge invariant 
(together with the Lagrangian). The presence of a ghost-like particle in the GS case renders the 
local description quite unappealing and for sure the best way to define the mechanism is just 
by adding the non-local counterterms. In the WZ case the local description is quite satisfactory 
and allows one to treat the bFF interaction as a real trilinear vertex, which takes an important 
role in the presence of a broken phase. The GS counterterms are, in practice, the same ones 
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(a) (b) (c) (d) 

Figure 6: The GS vertex with the non-local contributions for the BBB diagram symmetrized on each 
leg. 

as appearing in the analysis of axial QED with a non-local counterterm, as we are going to 
discuss next. 

3.2 Non-local count erterms: axial QED 

The use of non-local counterterms to cancel the anomaly is for sure a debated issue in quantum 
field theory since most of the results concerning the BRS analysis of these theories may not 
apply [H]. In the GS case we may ignore all the previous constructions and just require ab 
initio that the anomalous vertices are modified by the addition of a non-local counterterm that 
cancels the anomaly on the axial lines. 

Consider, for instance, the case of the BBB vertex of FigEJ where the regularization of the 
anomalies has been obtained by adding the three GS counterterms in a symmetric way [20J. 
In the BAA case only a single countertem is needed, but for the rest the discussion is quite 
similar to the BBB case, with just a few differences. These concern the distribution of the 
partial anomalies on the A and B lines in the case in which also BAA is treated symmetrically 
(equal partial anomalies). In this particular case we need to compensate the vertex with CS 
interactions, which are not, anyhow, observable if the A lines correspond to conserved gauge 
currents such as in QED. In this situation the Ward identities would force the CS counterterms 
to vanish. We will stick to the consistent definition of the anomaly in which only B carries the 
total anomaly a n and A is anomaly-free. The counterterm used in the GS mechanism both for 
BAA and BBB is nothing else but the opposite of the Dolgov-Zakharov (DZ) expression |21j, 
which in the BAA case takes the form 

C^ v (k, k h k 2 ) = -^ A e[/i, v, fci, k 2 \. (15) 

In the BBB case a similar expression is obtained by creating a Bose symmetric combination of 
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DZ poles, 

C AAA( k > fc l> fc 2) = - Q^-A: A e[/i, Z/, fci, fa] + e k A , fc 2, fc] + ^ fc 2 £ [^> ^> fc > fc l] j • ( 16 ) 

We have denoted by k the incoming momenta of the axial-vector vertex and by fa and fa the 
outgoing momenta of the vector vertices. We keep this notation also in the AAA case, since k 
will denote the momentum exchange in the s-channel when we glue together these amplitudes to 
obtain an amplitude of BIM type; this, we will analyze in the next sections. These expressions 
are consistent with the following equations of the anomaly for the BAA triangle 

kinC A y V (k, k±, fa) = 0, 
k 2v C A \^ v (k,fa,fa) = 0, 



(17) 



and for the BBB anomalous triangle 



k i» C AAA( k , k 2) = -y e [^> k, k 2 ], 
k 2vC% l A A (k,k 1 ,k 2 ) = -ye[A,/i, k, fa], 
k \ C AAA( k , k u h) = -ye[/^, v, fa, fa}. 

(18) 

So we can define a gauge invariant triangle amplitude, in both the BBB and BAA cases, 

by 

^AAA^{ k i k ^i k z) = ^AAa(^' ^1' ^2) + CaAa(^' ^1) ^2) 
^AVV S ( k ^ ^1, ^2) — ^AVvi k -> ^1' ^2) + C A yy{k, fa, fa). 

(19) 

Notice that the (fermionic) triangle diagrams, in the symmetric limit k\ — k\ — k\, is exactly 
the opposite of the DZ counterterms, as we will discuss in the next section, 

A A yy{k\ = k\ = k 2 ) = -C AVV A AAA (kj = k\ = k 2 ) = ~C AAA , (20) 

so the cancellation is identical at that point (only at that point), and the two vertices A GS 
vanish. It is rather obvious that the cancellations of these poles in BIM amplitudes, corrected 
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Figure 7: BIM-type amplitude with full GS vertices. For on-shell external lines the contributions 
from the extra poles disappear. 

by the GS counterterms, is identical only for on-shell external gauge lines. It is also quite 
straightforward to realize that the massless 5-model with the GS vertex correction is equivalent 
to axial QED corrected by a non-local term [5] that is described by the Lagrangian 



This theory is equivalent to the (local) formulation given in Eq. [5] and in [14J, where the 
transversality constraint (dB = 0) is directly imposed on the Lagrangian via a multiplier. 



Figure 8: Embedding of the BIM amplitude with GS vertices in a fermion/antifermion scattering 

Unitarity requires these DZ poles in the C counterterms to disappear from a physical am- 
plitude. To show that this is not the case, in general, consider the diagrams depicted in Figs. 
[7] and [SJ The structure of the GS vertex is, for BBB, given by (TTUT) with the three massless 
poles generated by the exchange of the pseudoscalar on the three legs, as shown in Fig. [6j For 
on-shell external lines, in this diagram the contributions from the extra poles cancel due to 
the transversality condition satisfied by the polarizators of the gauge bosons. However, once 
these amplitudes are embedded into more general amplitudes such as those shown in Fig. ( I8ll9lh 
the different virtualities of the momenta of the anomaly diagrams do not permit, in general, 
cancellation of the DZ extra poles introduced by the counterterm. 

We can summarize the issues that we have raised in the following points. 

1) The GS and the WZ mechanism have different formulations in terms of auxiliary fields. 





plus the counterterm 




(22) 
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Pi p'l 




(b) (c) 



Figure 9: The structure of the fermionic scattering amplitudes with spurious massless poles (graphs 
b) and c)). 

2) Previous analyses of axial gauge theories, though distinct in their Lagrangian formulation, 
are all equivalent to axial QED plus a non-local counterterm. The regularization of the gauge 
interactions, in these theories, coincides with that obtained by using the GS counterterm on 
the gauge vertex. In particular, the massless poles introduced by the regularization are not 
understood in the context of perturbative unitarity. 

A special comment is needed when we move to the analysis of Higgs- axion mixing. This has 
been shown to take place after electroweak symmetry breaking for a special class of potentials, 
which are not supersymmetric. The axion, which in the Stuckelberg phase is essentially a 
Goldstone mode, develops a physical component and this component appears as a physical 
pole. It is then clear, from the analysis presented above, that the regularization procedure 
introduced by the GS counterterm involves a virtual massless state and not a physical pole. 
This is at variance with the WZ mechanism, in which the bFF vertex is introduced from the 
beginning as a vertex and not just as a virtual state. In this second case, b can be decomposed 
in terms of a Goldstone mode and a physical pseudoscalar, called in [7j the axi-Higgs, which 
takes the role of a gauged Peccei-Quinn axion [7j. This is entitled to appear as a physical 
state (and a physical pole, massless or massive) in the spectrum. The re-formulation of the GS 
counterterm in terms of a pseudoscalar comes also at a cost, due to the presence of a ghost 
(phantom) particle in the spectrum, which is absent in the WZ case. There are some advantages 
though, since the theory has, apparently, a nice ultraviolet behavior, given the gauge invariance 
condition on the vertex. The presence of these spurious poles requires further investigations to 
see how they are really embedded into higher order diagrams and we hope to return to this point 
in the near future in a related work. There is another observation on this issue that is worth to 
mention: the anomaly is also responsible for a UV/IR conspiracy which is puzzling on several 
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grounds. For instance, the linear divergent terms A\ and A 2 in the Rosenberg representation 
[22J of the anomaly diagram are closely related to the infrared anomaly poles in the amplitudes 
A3 and Aq in the chiral limit, due to the Ward identity [23]. @ 

In the next section we will move to the analysis of the unitarity bound in the WZ case. 
As we have shown, in this case it is possible to characterize it explicitly. We will work in a 
specific model, but the implications of our analysis are general and may be used to constrain 
significantly entire classes of models containing WZ interactions at the LHC Before coming 
to the specific phenomenological applications we elaborate on the set of amplitudes that are 
instrumental in order to spot the bad high energy behavior of the chiral anomaly in s-channel 
processes: the BIM amplitudes. 

4 BIM amplitudes, unitarity and the resonance pole 

The uncontrolled growth of the cross section in the WZ case has to do with a certain class 
of amplitudes that have two anomalous (AVV or AAA ) vertices connected by an s-channel 
exchange as in Figj3]a). We are interested in the expressions of these amplitudes in the chiral 
limit, when all the fermions are massless. Processes such as AA — > A A, mediated by an anoma- 
lous gauge boson B, with on-shell external A lines and massless fermions, can be expressed in 
a simplified form, which is, also in this case, the DZ form. We therefore set k\ = k\ = and 
rrif = 0, which are the correct kinematical conditions to obtain DZ poles. We briefly elaborate 
on this point. 

We start from the Rosenberg form of the AVV amplitude, which is given by 

T XfJ,u = A 1 e[k u A, fi, u] + A 2 e[k 2 , A, fi, u) + A^e[k^ k 2 , u, A] 

+A 4 k^e[ki, k 2 , u, A] + A 5 k^e[ki, k 2 , fi, A] + A 6 k 2 e[ki, k 2 , /i, A] , (23) 

and imposing the Ward identities we obtain 

A x = k 2 2 A A + k x ■ k 2 A 3 
A 2 = k\A b + kx ■ k 2 A 6 
Mhufo) = -Ae(k u k 2 ) 

A 4 (k 1 ,k 2 ) = -A 5 {k 1 ,k 2 ), (24) 

4 We thank A. R. White for clarifying this point to us and for describing his forthcoming work on this issue. 
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where the invariant amplitudes A 3 , . . . , A 6 are free from singularities. In this specific kinematical 
limit we can use the following relations to simplify our amplitude 

2 

e[k 2 , A, //, v\ — - (k^elkt, k 2 , v, /i] + k^e^i, k 2l p, A] + k^efa, ^2, v, A]) 
2 

e\ki, A, /i, z/1 = — (fcMfci, k 2 , u, d + k\e\k^ k 2 , p, A] + fcMfti, fc 2 , z/, A]) (25) 
s 

where /c 2 = (k\ + A^) 2 = s is the center of mass energy. These combinations allow us to re- write 
the expression of the trilinear amplitude as 

T^ x = A 6 k x e[h, k 2 , v, /i] + (A 4 + A 6 ) (k'elh, k 2 , p, A] - k^h, k 2 , u, A]) . (26) 

It is not difficult to see that the second piece drops off for physical external on-shell A lines, 
and we see that only one invariant amplitude contributes to the result 

T^ A = A{(s)(ki + k 2 ) x s [h, k 2 , u, p] . (27) 

There are some observations to be made concerning this result. Notice that A 6 multiplies a 
longitudinal momentum exchange and, as discussed in the literature on the chiral anomaly 
in QCD [211 I2H [25], brings about a massless pole in s. We just recall that A 6 satisfies an 
unsubtracted dispersion relation in s at a fixed invariant mass of the two photons, (kf = k\= p 2 ) 

Aq(s,j)) = —\ dt (28) 

Jim) t - s 

and a sum rule 

f°° 1 

/ ImA 6 (t,p 2 )dt = — , (29) 

Jim} ^ 

while for on-shell external photons one can use the DZ relation [2T] 



ImA 6 (k 2 ,0) = -S(k 2 ), (30) 

7T 

to show that the only pole of the amplitude is actually at s = k 2 = 0. It can be simplified using 
the identity 

H 2 (l-a) + Li 2 (l-a- 1 ) = ~log 2 (a), (31) 

with 



a= Fj — Pf = \l-A^ } (32) 
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to give 

,2 



2vr 2 s 



Pf 



(33) 



We can use this amplitude to discuss both the breaking of unitarity and the cancellation of 
the resonance pole in this simple model. The first point has already been addressed in the 
previous section, where the computation of the diagrams in FigfS] has shown that only a BIM 
amplitude survives in the WZ case in the scattering process AA —>■ AA. If the sum of those 
two diagrams gives a gauge invariant result, with the exchange of the Z' described by a Proca 
propagator, there is a third contribution that should be added to this amplitude. This comes 
from the exchange of the physical axion x- We recall, in fact, that in the presence of Higgs- 
axion mixing, when b is the sum of a Goldstone mode and a physical axion x, each anomalous 
Z' is accompanied by the exchange of the x [I] • This is generated in the presence of electroweak 
symmetry breaking, having expanded the Higgs scalar <ft around a vacuum v 

<j) = -j= (v + X + t<f) 2 ) , (34) 

with the axion b expressed as linear combination of the rotated fields Xb an d G B 

o = aiXB + ot 2 G B = Xb + TT G b- (35) 

Mb Mb 

We also recall that the gauge fields get their masses M B through the combined Higgs- 
Stiickelberg mechanism, 



M B = yjM? + {q B g B vf. (36) 

In the phenomenological analysis presented in the next sections the contribution due to x nas 
been included. Therefore, in the WZ case, the total contributions coming from the several BIM 
amplitudes related to the additional anomalous neutral currents should be accompanied not 
only by the set of Goldstone bosons, to restore gauge invariance, but also by the exchange of 
the axi- Higgs. 

The cancellation of the resonance pole for s = M B is an important characteristic of BIM 
amplitudes, which does not occur in any other (anomaly-free) amplitude. This cancellation is 
the result of some amusingly trivial algebra which we reproduce just for the sake of clarity. 
Given a BIM amplitude and Proca exchange, we have 

Abjm = v, h, h] — ^_ O _ ^-k'W, v>, k[, k>] 



k 2 ' A ' ^ k 2 - Ml V J k 2 
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a n . r (71 * k X (Mf k 2 ) a n ,\'\r/ I ii ill 

= ^ £ i^ k ^ k ?\ k 2_ M 2 Mfi,v ,k,,k 2 } 

Q / (X 

= -jj-4v,v,Kk 2 ] jf^W,v', k i, k 2\- ( 37 ) 

This result implies that the amplitude is described - in the chiral limit and for massless external 
states - by a diagram with the exchange of a pseudoscalar (see Figj2]b)) and that the resonance 
pole has disappeared. It is clear, from fl37|) . that these amplitudes break unitarity and give a 
contribution to the cross section that grows quadratically in energy (~ s). Therefore, searching 
for BIM amplitudes at the LHC can be a way to uncover the anomalous behavior of extra 
neutral (or charged) gauge interactions. There is one thing that might tame this growth, and 
this is the exchange of the physical axion. We will show, working in a complete brane model, 
that the exchange of the x does lower the cross section, but insignificantly, independently from 
the mass of the axion. 



5 A realistic model with WZ counterterms 

Having clarified the relation between the WZ and GS mechanisms using our simple toy-model, 
we move towards the analysis of the issue of unitarity violation in the WZ Lagrangian at high 
energy. Details on the structure of the effective action of the complete model that we are going 
to analyze can be found in [U]. We just mention that this is characterized by a gauge structure 
of the form SU(3) x SU{2) x U(l) x U(1) B , where the U(1) B is anomalous. We work in the 
context of a two-Higgs doublet model with H u and Hd [Zj. In our analysis our setup is that of 
a complete model, in the sense that all the charge assignments are those of a realistic brane 
model with three extra anomalous U (1), but we will, for simplicity, assume that only the lowest 
mass eigenvalue taking part is significant, since the remaining two additional gauge bosons are 
heavy and, essentially, decoupled. 

We recall that the single anomalous gauge boson, B, that we consider in this analysis is 
characterized by a generator Yg, which is anomalous (TrFj 7^ 0) but at the same time has 
mixed anomalies with the remaining generators of SM and in particular with the hypercharge, 
Y. In the presence of a single anomalous U(l), here denoted as U(1)b, both the Z and the 
(extra) Z' gauge boson have an anomalous component, proportional to B. We also recall that 
the effective action of the anomalous theory is rendered gauge invariant using both CS and 
Green-Schwarz counterterms, while a given gauge invariant sector involves the exchange both 
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of the anomalous gauge boson and of the axion in the s-,t- and w-channels. 

In a previous work [TU] it has been shown how the trilinear vertices of the effective WZ 
Lagrangian can be determined consistently for a generic number of extra anomalous U(l). 
Here, the goal is to identify and quantify the contributions that cause a violation of unitarity 
in this Lagrangian. For phenomenological reasons it is then convenient to select those BIM 
amplitudes that have a better chance, at the experimental level, to be measured at the LHC 
and for this reason we will focus on the process gg — > 77. The gluon density grows at high 
energy especially at smaller (Bjorken variable) x-values. We choose to work with prompt final 
state photons for obvious reasons, the signal being particularly clean. To begin with, we will be 
needing the expressions of the Z77 and the Zgg vertices. In the presence of three anomalous 
£7(1), here denoted as U(1)b, both the Z and the (extra) Z' gauge boson have an anomalous 
component, which is proportional to the the anomalous gauge bosons of the interaction 
eigenstate basis (i=l,2,3). The photon vertex is given by [10] 



(zm)\ mf¥ c = -\zt*M £ [g Y e YYY R YYY + tiey ww R%Z w + 9 Y gle Y f ww R Y z ^ 

f 

+9Y92<7f ti Zal + / J 9B i gY92Vf R 2 



, „2 nBiYY TjBiYY . 2nBiWW ryBiWW 

+ 2_^9B t g Y Vf K Zny +gB t g 2 ff R; 



AfcK ± 0). 



(3* 



with I — 1,2,3 enumerating the extra anomalous neutral currents. The explicit expressions of 
the rotation matrix A can be found in [9j [10]. We have defined 

R YYY = (0 A ) YZl (0 A f Yv R™ w = (0 A ) W[iZl (0 A f WsV (39) 

and the triangle A^yy(m./ 7^ 0) is given by [10] 

A Avv( m f 1 0, k u k *) = ^i dx i d V^) 
{e[ki, A, fi, u] [y(y - \)k\ - xyk x ■ k 2 ] 
+e[k 2 , A, /i, v\ [x(l - x)k\ + xyki ■ k 2 ] 
+e[ki, k 2 , A, v\ [x(x — l)k^ — xyk%] 
+e[k 1} k 2 , A,/i] [xyk\ + (1 - y)yk v 2 ]} , 

A(m f ) =m 2 f + x(x - l)k\ + y(y - \)k\ - 2xyh ■ k 2 . (40) 
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We have defined the following chiral asymmetries 




- Qkf( T l f r, 

= (Tiff, 
Q' ! J Qy./lh. , 



(41) 



with Q B and Q Y denoting the charges of the chiral fermions and T| is the generator of the 
third component of the weak isospin, while the R factors are products of A matrix elements. 
The matrix A relates the interaction eigenstate basis of the generators (Y B , Y, T 3 ) to those of 
the mass eigenstate basis (T z , Tz>, Q), of the physical gauge bosons of the neutral sector, Z, Z' 
and Ay. They are given by 



R Y z^ = H(o A ) YZl (o A r Yl ] 

R Y z^ w = [2(o A ) W31 (o A ) YZl (o A ) Yl + (o A y w jo A ) YZl ] 

R^=[3(o A ) BiZl (o A y W31 ] 

R Y ™ = [2(O a ) YZi (O a ) Yi (O a ) W31 + [0 A ) W3Zl {0 A f Yl ] 

Rl zni = (° A )yj( C ' A )b 1 Zi 

R Z ^ W = [(0 A )w 31 {0 A ) BiZl ] 

Rz^ W = [2(O a ) BiZi (O a ) W31 (O a ) Yi ] . 



These expressions will be used extensively in the next section and computed numerically in a 
complete brane model. 

5.0.1 MLSOM with one Higgs doublet and an extra singlet 

Another possible framework for the Higgs sector is to consider a Lagrangian containing one 
SU(2) W Higgs doublet H u and an extra singlet 0. 

In this context all the features concerning the mass matrix for the gauge bosons remain the 
same, in fact the covariant derivatives act on the fields as follows 



(42) 




(43) 
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where and are respectively the charges of the higgs fields under the extra anomalous 
U(l) symmetry. 

After the spontaneous symmetry breaking, expanding around the vacuum of the two Higgs 
bosons, we have 

K = v. + BeH °+i hnH ", 4, = V , + M±pl. (44) 

By this procedure we obtain a 3 x 3 mass matrix in the mixing of the neutral gauge bosons, 
which is similar to that obtained in [9]. 

In the Higgs sector the structure of the Peccei-Quinn potential is similar to that obtained 
in the presence of two Higgs doublets. In fact, the symmetric potential is given by 

V PQ = fi 2 u HtH u + X uu (HlH u ) 2 - 2\ u ^HlH u )(t>*(j) + + A^(0» 2 , (45) 

where the coefficients ji u and have mass dimension 1 and X uu , A U( ^, are dimensionless, 
while the PQ breaking terms can be written as 

= Ji{HiH u )<f> e iq *^ + A(0*0)0 e ig "^ + c.c. , (46) 

where Jjl and A have mass dimension 1. The CP-odd sector is still characterized by a 3 x 3 
matrix similar to O x , which we call for simplicity 0' x ; this allows for the rotation from the 
interaction eigenstates basis to the physical basis, as follows 



/ Im#° 






( 


X 


\ 


Im0 




= 0' x 




G? 




V b 


J 




\ 


G° 


/ 



Therefore, in this picture the spectrum in this sector is unchanged, we have a physical axion 
and its mass exhibits a combination of two effects, the Higgs vevs and the presence of a PQ 
breaking potential. The parameters can be tuned to be small enough to have a light physical 
axion. In general, bounds coming from experimental, astrophysical and cosmological sources 
give an upper limit for the axion mass, which is around 10 meV j26j.[27j. 

5.1 Prompt photons, the Landau- Yang theorem and the anomaly 

Since the analysis of the unitarity bound will involve the study of amplitudes with direct 
photons in the final state mediated by a Z or a Z' in the s-channel, we briefly recall some facts 
concerning the structure of these amplitude and in particular the Landau- Yang theorem. The 
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Figure 10: Vanishing and non-vanishing amplitudes mediated by an axial-vector spin 1. 

condition of transversality of the Z' boson (ez> ■ k = 0) is essential for the vanishing of this 
amplitude. A direct proof of the vanishing of the on-shell vertex can be found in [TT|, [9] . 

The theorem states that a spin 1 particle cannot decay into two on-shell spin 1 photons 
because of Bose symmetry and angular momentum conservation. Angular momentum con- 
servation tells us that the two photons must be in a spin 1 state (which forces their angular 
momentum wave function to be antisymmetric), while their spatial part is symmetric. The 
total wave function is therefore antisymmetric and violates the requirement of Bose statistics. 
For these reasons the amplitude has to vanish. For a virtual exchange mediated by a Z' the 
contribution is vanishing -after summing over the fermions in each generation of the SM-, the 
theory being anomaly-free, in the chiral limit. The amplitude is no n- vanishing only in the 
presence of chiral symmetry breaking terms (fermion masses), which can be induced both by 
the QCD vacuum and by the Yukawa couplings of SM in the presence of electroweak symmetry 
breaking. For this reason it is strongly suppressed also in the SM. 

However, the situation in the case of an anomalous vertex is more subtle. The BIM am- 
plitude is non-vanishing, but at the same time, as we have explained, is non-resonant, which 
means that the particle pole due to the Z' has disappeared. For the rest it will break unitarity 
at a certain stage. 

In fact , a cursory look at the AVV vertex shows that if the external photons are on-shell and 
transverse, the amplitude mediated by this diagram is proportional to the momentum of the 
virtual Z , k^. This longitudinal momentum exchange does not set the amplitude to unless 
the production mechanism is also anomalous. We will show first that in the SM these processes 
are naturally suppressed, though not identically 0, since they are proportional to the fermion 
masses, due to anomaly cancellation. We start our analysis by going back to the AVV diagram, 
which summarizes the kinematical behavior of the Z77 amplitude. 

Let k\ and &2 denote the momenta of the two final state photons. We contract the AVV 
diagram with the polarization vectors of the photons, £1^ and e 2 i/ of the Z boson, e\, obtaining 
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a 1 f l r 1 ^ x X y 

exe^u^Avvih, k 2 ,m f ^ 0) = » dx dy— — - — 

Jo Jo rnj - 2xyh ■ k 2 

e\S lfi e 2 v |£A^a (A£ - ^l) + K k 2 (^aA/rf - S^^K)] , ( 47 ) 

where we have used the conditions 

Jv-^ — — fv<2 — — 

e^kl = e 2v k v 2 = 0. (48) 

It is important to observe that if we apply Schouten's identity we can reduce this expression 
to the form 

A^h, k 2 ,m f ^ 0) = e A (k x + k x ) T s {k\m s ) (49) 



with 



Tj{k 2 ,m s ) = J f (k 2 )e[k u k 2 ,Ei,e 2 ] 



(50) 



which vanishes only if we impose the transversality condition on the polarization vector of the 
Z boson, e\k x = 0. Alternatively, this amplitude vanishes if the anomalous Z-photon-photon 
vertex is contracted with another gauge invariant vertex, as discussed above. The amplitude 
has an anomalous behavior. In fact, contracting with the k\ four- vector we obtain 

a ( 1 f 1 f 1 ~ x k\k x xy \ 

kxe^e^A/Jyikuk^rrif ^ 0) = I -— dx dy ■ - ^—j- - — J e[k u k 2 ,e u e 2 ] 

h I dX I dy m 2 2 -2xvk^ko ) £ ^'^,£i,£ 2 ] 

+ — / dx dy— — — I e[fci, fc 2 ,ei,e 2 ] (51) 



7T^ 7T Z 







o mj — 2xyk 1 ■ k 2 



where the mass-independent and mass-dependent contributions have been separated. Summing 
over an anomaly-free generation, the first of these two contributions cancel. 
For this reason, it is also convenient to isolate the following quantity 

g f (k ,m f ) = - 1 dxj^ dy m) _ 2xykih e[ kl ,k 2 ^e 2 ] (52) 
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which is the only contribution to the triangle amplitude in the SM for a given fermion flavor /. 

We illustrate in Fig{TU] the contributions to the qq annihiliation channel at all orders (first 
graph), at leading order (second graph), and the BIM amplitude (third graph), all attached 
to an AVV final state. In an anomaly-free theory all these processes are only sensitive to 
the difference in masses among the flavors, since degeneracy in the fermion mass sets these 
contributions to zero. Only the second graph vanishes identically due to the Ward identity 
on the qq channel also away from the chiral limit. For instance in the third diagram chiral 
symmetry breaking is sufficient to induce violations of the Ward identity on the initial state 
vertex, due to the different quark masses within a given fermion generation. 

To illustrate, in more detail, how chiral symmetry breaking can induce the exchange of a 
scalar component in the process with two prompt photons, we start from the SM case, where 
the Z*77 vertex, multiplied by external physical polarizations for the two photons becomes 



V x 



92 



2cos# 



ir 



« 2 E(Q/)M//(fcV/) 
/ 



(53) 



and consider the quark antiquark annihilation channel qq — > Z* — > 77. We work in the parton 
model with massless light quarks. We can rewrite the amplitude as 



M 



ffz* lL \x,t; v z* lV 



where we have introduced the Zqq vertex 

Vffz = v( Pl )T x 'u(p 2 ) T A = ig z i x {g v - 9a1 5 ), 
and the expression of the propagator of the Z in the gauge 

k x k x ' 



n 



A,A',£ 



A' 2 



M\ l 



,AA' 



9 



k 2 - £MJ 



To move to the unitary gauge, we split the propagator of the Z as 



n 



—% 



A,A',£ 



k 2 



9 



XL* 



k x k 



+ 



—1 



k x k x ' 



Ml 



k 2 - £Mf 

and go to the unitary gauge by choosing £ — >• 00. The amplitude will then be written as 

k x k x " 



(54) 



(55) 



(56) 



(57) 



M = w(pi)r A M (p 2 



k 2 



Ml 



,AA 



Ml 



9 
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Ml 



k 



92 

2 cos 9 



e 2 J2(Qf) 2 9Af^(k 2 ,m f ) 



2 cos 9w 



f 



(58) 
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Figure 11: Anomalous contributions in gg — ► 77 



Clearly, at the Born level, using the Ward identity on the left Vffz* vertex, we find that 
the amplitude is zero. This result remains unchanged if we include higher order corrections 
(strong/electroweak), since the structure of this vertex is just modified by a Pauli (weak-electric) 
form factor and the additional contribution vanishes after contraction with the momentum k x . 
This amplitude is however non- vanishing if we replace the Vf/z* vertex with a V gg z* vertex, 
where now we assume that the new vertex is computed for non-zero fermion masses (i.e. away 
from the chiral limit). In this case we use the Ward identity 



k V pvp 



(Pl+P2)pG 



pufi 



z 92 



2 cos 9 W 



2 vp,af3 



Pl*P2f3 



7T 2 7T 2 jq 



dxi 



1 — x\ 



dxo 



(59) 



with 

A q = -xyk 2 + m 2 q) (60) 

where the constant term (1/V 2 ) vanishes in an anomaly- free theory. It is convenient to 
define the function 



G q (k 2 ,m q ) = e[e 1 g,€2g,Pi,P2]m 2 q l q (m q ) 



(61) 



with 



l g (m q ) 



-4 / dxi J dx 2 



71 JO 



(62) 



where e is the polarization of the gluon, which allows one to express the squared amplitude as 



e 2 92 



2cos9 w J 4M| 



E 

1 



9 q AQ 2 Mm q )) \Y,£<%Gf{mj) \ , 



(63) 
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with s = k 2 . Notice that in the large energy limit 1 q ~ Jj ~ 1/s 2 [22] • This shows that double 
prompt photon production, in the SM, is non-resonant and is proportional to the quark masses, 
neglecting the contributions coming from the masses of the leptons. 



6 Gauge parameter dependence in the physical basis 



When this analysis is extended to a complete anomalous model such as the MLSOM [TJ [9], [TO] 
even the direct proof of the cancellation of the gauge dependence in the Z' exchange is quite 
complex and not obvious, although it is expected at a formal level. We recall that we are 
working with a broken phase and the axion has been decomposed into its physical component 
(x, which is the axi-Higgs) and the Goldstone modes of the extra Z', Gz 1 - In pjj] it has been 
shown that the counterterms of the theory can be fixed in the Stuckelberg phase and then 
re-expressed, in the Higgs-Stuckelberg phase, in the physical base. Although the procedure 
is formally correct, the explicit check of cancellation of these gauge dependences is far from 
trivial and is based on some identities that we have been able to derive after some efforts, which 
confirm the correctness of the approach followed in 0, [10] for the determination of the effective 
Lagrangian of the model after electroweak symmetry breaking. 

The matrix O x , needed to rotate into the mass eigenstates of the CP-odd sector, relating 
the axion x an d the two neutral Goldstones of this sector to the Stuckelberg field b and the 
CP-odd phases of the two Higgs doublets satisfies the following relation 



/ ImEl 


\ 




f 


X 


\ 


ImH° d 




= O x 




G? 




\ h 


) 




\ 


Gl 


J 



where the Goldstones in the physical basis are obtained by the following combination 
G z = G? 
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Here we have defined the following coefficients 

fu = g 2 0^ W3 - g Y 0^ Y - q*g B Oi B f d = g 2 0$ W3 - g Y 0^ Y - q B d g B 0^ B 
fu,B = g20z>w 3 - gy°z>Y - Qu gBOz>B fd.B = g20z'w 3 - gvOziy - q B g B Oj, B , (65) 

and the q^ d charges are defined in Table ([2]). The relations containing the physical Goldstones 
can be inverted so we obtain 

G° = dG z + C 2 G Z ' 

G° 2 = C[G Z + C' 2 G Z ' , (66) 

where we give the explicit expression only for the coefficient C[, since this is the one relevant 
for our purposes. Then, after the orthonormalization procedure, we obtain 

C[ = ; ° 2 . (67) 

V4 + c 'i 

We illustrate the proof of gauge independence from the Z gauge parameter of the amplitudes 
in Figfl2l In this case the cancellation of the spurious poles takes place via the combined 
exchange of the Z propagator and of the corresponding Goldstone mode Gz- If we isolate the 
gauge-dependent part in the Z boson propagator we obtain 



+4 x (Ac^e\p 1 ,p 2 ,n,u}) k2 _ (4c^e[k 1 ,k 2 ,fJ,',v'}) 

+4 x (Ac z 7 e[ Pl ,p 2 ,fi,u]) k2 _ ^e^,^, //>']) , (68) 
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where 



77 



(69) 



The coefficient c^ 7 in the M. b% amplitude (Gz exchange) must be compared with the massless 
coefficients c{ of the A4a% amplitude (Z boson exchange) and the explicit expressions of the 
coefficients Cyy and F are worked out in the next section. 
Adding the contributions of the two diagrams we obtain 



Mas + M B a 



n - [g B g Y e? YY RZ Y + 9B9 2 2 9f ww RZ w ] e[ Pl ,p 2 , u]} ^~ 

{ y [9B 9Y ef YY R B z YY + 9B9lef ww RZ w ] **> ^ A } 



+4 x 
4 



W-tzMl 



(70) 



At this point, the pattern of cancellation can be separated in three different sectors, a pure 
BYY sector, a pure BWW and mixed BYY-BWW sectors, and it requires the validity of the 
relations 



(%9B&er~RgZ) 2 ± + 4 {^{0 Yl fO^ = 0, (71) 
(I « 2 2 rC) 2 ^1 + 4 (^O^fO^ 2 = 0, (72) 

(angs^ef^R^^gsdef^R^) ± + 4 (^(^ 7 ) 2 ^ «) 2 3 *^ = 0. 

(73) 

We have been able to verify that these relations are automatically satisfied because of the 
following identity, which connects the rotation matrix of the interaction to the mass eigenstates 
A to a component of the matrix O x . This matrix appears in the rotation from the basis of 
Stiickelberg axions to the basis of the Goldstones Gz and Gz< and of the axi-Higgs x- The 
relation is 

0^ Z -L = 20^C[^-, (74) 
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with 
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(75) 
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with Mi the Stiickelberg mass. The origin of this connection has to be found in the Yukawa 
sector and the condition of gauge invariance of the Yukawa couplings. 

7 Unitarity bounds: the partonic contribution gg —> 

In this section we perform the analytical computation of the cross section for the process 
99 ~~ * 77 with the two on-shell gluons (g) and the two on-shell photons in the final state. 
The same computation is carried out both in the SM and in the MLSOM of [7J, where the 
charge assignments have been determined as in |12j . to determine the different behavior of 
these amplitudes in the two cases. The list of contributions that we have included are all shown 
in FigJTTl We report some of the results of the graphs in order to clarify the notation. For 
instance we obtain for diagram (a) 




(76) 



where we have defined the coefficients 



A = 9l [9re Y q 0^ z + 920? 0* z + QB^Oiz] , 




(77) 



and the mass of the extra Z' is expressed as 




/ B2 2 , B2 2\ 2 I B 2 , B 2\ 

[Qu v u + q d v d )g B , x B = (q u v u + q d v d )g B . 
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We have also defined v = \Jv 2 + v%, where v u and v a are the vevs of the two Higgs bosons, in 
the scalar potential J7] and g 2 = ' g\ + g\. We recall that A((s) is approximately given by 



1 



2w 2 s 2s 2 



m 2 t f m 2 f ( m 2 t 
-4 + [-L log [-L 



(80) 



at large values of s. We have seen that the SM contribution, in the presence of a massive fermion 
circulating in the loop is suppressed by a factor that is O (m 2 / s 2 ) . In the case of an anomalous 
model this contribution becomes subleading, the dominant one coming from the anomalous 
parts, proportional to the chiral asymmetries 6f of the anomalous charge assignments between 
left-handed and right-handed fermion modes. The amplitude in diagram (b) is given by 
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where, for convenience, we have defined 
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We have defined the model-dependent parameters [9] 
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Proceeding in a similar way, graph c) gives 
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where we have used Eq.( l3T|) and we have set 
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Using the Yukawa couplings shown in [9J it is convenient to define the following coefficients, 
which will be used below 

c gg = 9lc x ' g } q = u,d 

c x ' y f = e 2 c x ' f , q = u,d, u, e. (86) 

We have used a condensed notation for the flavors, with u = {u, c, t}, d = {d, s, b}, v = {z/ e , 
v^, v T } and e = { e, fi, r}. The couplings of the physical axion to the fermions are given by 
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We have also relied on the definitions of O x introduced in a previous work [9] 



and the fermion masses have been expressed in terms of the Yukawa couplings by the relations 

v u T u v d T d v u T u v d T e 

m u = — — , m d = m v = — ;=-, m e = (89) 



The cross section for the amplitude (d) in Fig. ffTTl) is given by 
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With these notations, we are now ready to express the cross section for graph e) as 
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Finally, the cross section for the Z' exchange is given by 
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Figure 13: Partonic cross section for the anomalous process gg — ► 77 with rrif = 0, tan/3 = 40, 
gs = 0.1, M x = 10 GeV and Mi = 800 GeV. The solid line refers only to the exchange of the Z and 
the Z', while the dashed line refers to the complete cross section including the x exchange. 

where, again, in order to simplify the notation we have denned the coefficients 
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In the SM case we work in the unitary gauge, being a tree-level Z exchange, and we have 
only one contribution whose explicit expression is given by 
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where we have denned the SM coefficients 
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Figure 14: The same as in Fig. [13] but with rrif / 0. 
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(95) 

The partonic SM cross section is in agreement with unitarity in its asymptotic behavior, which 
is given by 
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In our anomalous model the complete expression of the same cross section is given by 
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Figure 15: Partonic cross section for the anomalous process gg — ► 77, with tan/3 = 40, gs = 0.1, 
M x = 10 GeV and Mi = 800 GeV. The solid line refers to the SM with the exchange of the Z boson. 
The dashed line refers to the MLSOM case. The point of minimum divides the anomaly-free region 
from the region where the anomalous contributions dominate. 
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Figure 16: Partonic cross section for the anomalous process gg — ► 77 with gs = 0.1, M x = 10 GeV 
and Mi = 800 GeV. The lines refer to the cross section evaluated for different values of tan (3. 
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At high energy we can neglect the mass of the axion (s — M 2 xs s) and from the limits 

A(s — * 00, rrif) = 1, Co(s — > 00, mj) = 0, (99) 
we obtain that the total cross section reduces to 

a MLSOM {s ^ Qo) = K?^ (10Q) 
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Figure 17: Partonic cross section for the anomalous process gg — > 77 with tan (3 = 40, M x = 10 GeV 
and Mi = 800 GeV. The lines refer to the cross section evaluated for different values of the coupling 
constant gs- 



which is linearly divergent and has a unitarity bound. /C = /C(s&, gs, cts(s), tan f3) is defined by 
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The derivation of the unitarity bound for this cross section is based, in analogy with Fermi 
theory, on the partial wave expansion 
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Figure 18: Partonic cross section for the anomalous process gg — > 77, with tan/3 = 40, gs = 0.1 and 
M x = 10 GeV. The lines refer to the cross section evaluated for different values of the Stiickelberg 
mass M\. 



with an s-wave contribution given by 
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Since unitarity requires that |/;| < 1 we obtain the bound 
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The bound is computed numerically by looking for values s& at which 
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where in the total parametric dependence of the factor /C, fC(s b ,gB, as(sb),tan/3), we have 
included the whole energy dependence, including that coming from running of the coupling (up 
to three- loop level). We will analyze below the bound numerically in the context of the specific 
brane model of [12], |28j . 

8 Couplings and Parameters in the Madrid model 

We now turn to a brief illustration of the specific charge assignments of the class of models that 
we have implemented in our numerical analysis. These are defined by a set of free parameters, 
which can be useful in order to discern between different scenarios. In our implementation 
we rotate the fields from the D-brane basis to the hypercharge basis and at the same time we 
redefine the abelian charges and couplings. The four U(l) in the hypercharge basis are denoted 
U(l)xi with i = A,B,C and Uy, where the last is the hypercharge £7(1), which is demanded 
to be anomaly-free. This fixes the hypercharge generator in the hypercharge basis in terms of 
the generator q a (a = a, b, c, d) in the D-brane basis. The U(l) a and U(l)d symmetries can be 
identified with (three times) the baryon number and (minus) the lepton number respectively. 
The U(l) c symmetry can be identified with the third component of right-handed weak isospin 
and finally the £7(1)& is a PQ-like symmetry. Specifically the hypercharge generator is given by 

Y = ^(qa + 3q d )-\qc, (108) 

which in fact is a linear combination of the two anomaly-free generators (q a + 3qd) and q c , while 
the orthogonal combinations 

X A = 3g a - q d , X B = q b , (109) 

represent anomalous generators in the hypercharge basis. Note that relation (11081) must be 
imposed in these models in order to obtain a correct massless hypercharge generator as in 
the SM. The set (Y, A, B) does not depend on the model, while the fourth generator X c is 
model-dependent and is given by 

/ 3 @2 n a2 3 p j3 2 \ 
X c = I — ^ q a + 6pn bl q b + 2n cl q c H q d I . (110) 

As can be seen in the detailed analysis performed in [T2] [28], the general solutions are parametrized 
by a phase e = ±1, the Neveu-Schwarz background on the first two tori $ = 1 — 6, = 1, 1/2, 
the four integers n a 2,n b \,n c \ and n^ 2 , which are the wrapping numbers of the branes around 
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the extra (toroidal) manifolds of the compactification, and a parameter p = 1,1/3, with an 
additional constraint in order to obtain the correct massless hypercharge 

n a = 777r( n a2 + 3pn d2 ). (Ill) 

This choice of the parameters identifies a particular class of models which are called Class A 
models [28J and all the parameters are listed in Tables ffl~ir2f3]) . Whether anomalous or not, the 
abelian fields have mass terms induced by the Stiickelberg mechanism, the mass matrix of the 
U{1) gauge bosons in the D-brane basis is given by the following expression 

3 

(M 2 U = g a gpM 2 s J2 c ?^i ( 112 ) 
i=i 

where Ms is some string scale to be tuned. Greek indices run over the D-brane basis {a, b, c, d}, 
the Latin index % runs over the three additional abelian gauge groups, while the g a and gp are 
the couplings of the four U(l). The eigenvectors u>j (i=Y,A,B,C) and their eigenvalues Aj for 
the matrix (M 2 ) Q/ g have been computed in terms of the various classes of models in reference 

m 

Wy = (113) 



Wi = t^t {w ia ,w ib ,w ic , 1} (114) 

where the Wi = A,B,c are the components of the eigenvectors. On the basis of the analysis of 
the mass matrix one can derive a plot of the lightest eigenvalue M.3, which corresponds to 
the X B generator in the hypercharge basis [2B]- We have reproduced independently the result 
for the lightest eigenvalue M 3 , which is in good agreement with the predictions of [2B]- We 
have implemented numerically the diagonalization of the mass matrix and shown in Fig{19] the 
behavior of the ratio M3/M5 as a function of the wrapping number n a2 and for several values of 
the ratio 1Z = gd/g c - This ratio, which characterizes the couplings of U(l) c and U(l)d, appears 
as a free parameter in the gauge boson mass matrix. Ms, the string scale, is arbitrary and can 
be tuned at low values in the region of a few TeV. 

We have selected a Stiickelberg mass M 3 (which is essentially the mass of the extra Z') of 
800 GeV. From Figfl9l if we choose the curve with 1Z = 1 at the peak value, then this values 
is 13% of the string scale, which in this case is lowered to approximately 6.1 TeV. It is quite 
obvious that the mass of the extra Z' can be reasonably assumed to be a free parameter for all 
practical purposes. 
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Figure 19: Lightest eigenvalue related to the generator Xb for a model of Class A. We have chosen 
the couplings g c and g d such that their ratio 1Z = gd/g c - 
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Table 1: Parameters for a Class A model with a D6-brane . 



The matrix Ei a = (wi) a constructed with the eigenvectors of the mass matrix defines the 
rotation matrix Ui a for the U(l) charges from the D-brane basis {a, b, c, d} into the hypercharge 
basis {Y, A, B, C} as follows 

^ffa, U ia = 9 —E iai (i = Y,A,B,C). (115) 
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So for the hypercharge we find that 
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to be identified with the correct hypercharge assignment given in expression (11081) . This iden- 
tification gives a relation between the gauge couplings 
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Table 2: Higgs charges in the Madrid model. 
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Table 3: SM spectrum charges in the .D-brane basis for the Madrid model. 



For the third generator a similar argument gives 

9a i 9b i 9c 9d / 1 1 o \ 

Qb = — w Ba Qa + —w Bb q b + — w Bc q c + — w Bd q d = q b , (118) 
9b 9b 9b 9b 

where we identify the gauge symmetry B corresponding to the lightest mass eigenvalue as an 
anomalous generator q B - The charges q B of the SM spectrum are given in Table HI We recall 
that given a particular non-abelian SU (N) gauge group, with coupling g N , arising from a stack 
of N parallel branes and the corresponding U(l) field living in the same stack, the two coupling 
constants are related by g\ = g^/V^N. Therefore, in particular the couplings g a and g b are 
determined using the SM values of the couplings of the non-abelian gauge groups 

91 = ^, 9l = g \, (119) 

where g^ is the SU{2) L gauge coupling, while g c and are constrained by relation (j!17p . 
Imposing gauge invariance for the Yukawa couplings [0] we obtain the assignments for the 
Higgs doublets shown in Table (j2J). 



8.1 Anomalous and anomaly-free regions: numerical results 

We have implemented in Candia [29j a numerical program that will provide full support for the 
experimental collaborations for their analysis of the main signals in this class of models at the 
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Table 4: Fermion spectrum charges in the Y-basis for the Madrid model [12]. 
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Figure 20: Partonic cross section for the anomalous process gg — > 77 with tan/? = 40, gs = 0.1, 
M x = 10 GeV and Mi = 1 TeV. The lines refer to the cross section evaluated for different values of 
the center of mass energy y/s. 

LHC, with the charge assignments discussed in the previous section. Candia has been planned 
to deal with the analysis of extra neutral interactions at hadron colliders in specific channels, 
such as Drell-Yan processes and double prompt photon processes with the highest precision, 
and it is under intense development. The program is entirely based on the theory developed in 
[3, EH El E] and in this work and is tailored to determine the basic processes, which provide the 
signal for the anomalous and anomaly-free extra Z' both from string and GUT models. The 
QCD corrections, including the parton evolution, is treated with extreme accuracy using the 
theory of the logarithmic expansions, developed in the last several years. Here we will provide 
only parton level results for the anomalous processes that we have discussed in the previous 
sections, which clarify the role played by the WZ Lagrangian in the restoration of unitarity at 
high energy. A complete analysis for the LHC is under way. 
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Figure 21: The partonic cross section with a parameter choice as in Fig. [20l The dashed upper lines 
refer to MLSOM cross sections evaluated with and without the x exchange, while the decreasing solid 
line refers to the SM. 



• a reduction by the exchange of \ 

We show in Fig. [13] a plot of the (small) but increasing partonic cross section for double 
prompt photon production from anomalous gluon fusion. We have chosen a typical SM- 
like value for the coupling constant of the extra Z' included in the analysis and varied 
the center of mass energy of few GeV's around 4.2 TeV. We show two plots, both in the 
brane model, one with the inclusion of the axi-Higgs x an d one without it, with only 
the exchange of the Z and Z'. Notice that the exchange of the x is a separate gauge 
invariant contribution. We have chosen a Stiickeiberg mass of 800 GeV. The plots show 
the theoretically expected reduction of the linear growth of the cross section, but the 
improvement is tiny, for these values of the external parameters. In these two plots, the 
fermion masses have been removed, as we worked in the chiral limit. The inclusion of 
all the mass effects in the amplitude has an irrelevant effect on the growth of the cross 
section. This is shown in Fig. [14] where, again, the inclusion of the axi-Higgs lowers the 
growth, but only insignificantly. We have analyzed the behavior of the cross sections in 
the presence of a light (about ten or more GeV's) axi-Higgs, but also in this case the 
effects are negligible. This feature can be easily checked from Eq. [97] in fact, the mass 
term M x is contained in the denominator of the propagator for the scalar, and in the 
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Figure 22: Partonic cross section plotted for different values of the coupling constant gs- The 
parameters are chosen as before. 

TeV's region we have (s — M?) ~ s. The numerical value of the unitarity bound remains 
essentially unchanged. 

• Anomaly-free and anomalous regions 

An interesting behavior shows up in Fig. [T5J where we compare the results in the SM 
and in the MLSOM for the same cross section, starting at a lower energy. It is clear that, 
while the SM result, being anomaly-free, is characterized by a fast-falling cross section, in 
the MLSOM it is very different. In particular, one finds a region of lower energy, where 
essentially the model follows the SM behavior (below 1 TeV) - but smaller by a factor 
of 10 - the growth of the anomalous contributions still being not large; and a region of 
higher energy, where the anomalous contributions take over (at about 2 TeV) and which 
drive the growth of the cross section, as in the previous two plots. There is a minimum 
at about 1.2 TeV, which is the point at which the anomalous subcomponent becomes 
sizeable. 

• tan/3, gb and M\ variations 

The variation of the same behavior shown in the previous plot with tan/3 is shown in 
FigHH where we have varied this parameter from small to larger values (10-50). The 
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Figure 23: Partonic cross section for the anomalous process gg — > 77, tan/3 = 40, gs = 0.1, and 
M x = 10 GeV plotted for different values of the Stiickelberg mass M\. 

depth of the minimum increases as we increase this value. At the same time, the cross 
section tends to fall much steeper, starting from larger values in the anomaly-free region. 

In FigHT] we study the variation of a as we tune the coupling of the anomalous U(l), 
gs- A very small value of the coupling tends to erase the anomalous behavior, rendering 
the anomalous components subleading. The cross section then falls quite fast before 
increasing reaching the bound. The falling region is quite visible for the two values of 
gs = 0.001 and 0.01, showing that the set of minimum points, or the anomaly-free region, 
is pushed up to several TeV's, in this case above 4.5 TeV. The unitarity bound is weaker, 
being pushed up significantly. The situation is opposite for stronger values of gs- 

A similar study is performed in FigJT8l but for a varying Stiickelberg mass. As this mass 
parameters increase, the anomaly-free region tends to grow wider and the cross section 
stabilizes. For instance, for a value of the Stiickelberg mass around 5 TeV, the region 
in which a has a normal behavior moves up to 3.5 TeV. The explanation of this result 
has to be found in the fact that the anomalous growth is controlled by the mass of the 
anomalous Z' in the s-channel, appearing in the denominator of the cross section. This 
suppression is seen both in the direct diagram and in the counterterm diagram, which 
describes the exchange of the axi-Higgs. Obviously, it is expected that as we reduce 
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Figure 24: Partonic cross section for the anomalous process gg — > 77, tan/3 = 40 Mi = 800 GeV 
and M x = 10 GeV. The different plots show a comparison between the SM and the MLSOM cross 
sections, at very high energies, for small value of the coupling constant gs- 

the coupling of the anomalous gauge boson, the anomalous behavior is reduced as well. 
The different behavior of the cross section in the SM and MLSOM cases can easily be 
inferred from FigEHl having chosen a Stiickelberg mass of the order of 1 TeV. A similar 
behavior is quite evident also from Figj2Tl from which it appears that in the MLSOM the 
deviations compared to the SM partonic predictions get sizeable at parton level already 
at an energy of 4-6 TeV. Notice also from Fig HH that the presence of the axi-Higgs seems 
to be irrelevant for the chosen values of the couplings and parameters of the model. We 
show in Fig. [22] a plot of the dependence of the predictions on g% at larger energy values, 
which appears to be quite significant. 

Furthermore, in the TeV's region the MLSOM predictions for small values of the cou- 
pling constant (g B = 0.001) go below the SM prediction and this is due to the axi-higgs 
exchange, which is negative in this kinematical domain. Moving below to 1 TeV the 
axi-higgs interference has an opposite sign and the MLSOM predictions are above the 
SM. A similar analysis, this time for a varying Stiickelberg mass, is shown in Figf23l and 
also in this case, as in the previous one, the results confirm that this dependence is very 
relevant. 

Finally, in FigJ2H we plot the SM and MLSOM results on a larger interval, from 10 to 
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Figure 25: Behavior of the mass of the additional anomalous Z' as a function of tan (3 for different 
values of qb- The variations are very small. 



140 TeV, from which the drastically different behavior of the two cross sections are quite 
clear. Notice that as we lower gs, for instance down to 10~ 4 , the anomaly-free region 
extends up to energy values that are of the order of 200 TeV or so. We conclude that 
the enhancement of the anomalous contributions with respect to the SM prediction are in 
general quite large and very sensitive to the mass of the extra Z' and to the strength of 
the anomalous coupling. Interestingly, a very weakly coupled Z' gives a cross section that 
has a faster fall-off compared to the SM case in the anomaly-free region. The MLSOM 
and SM predictions intersect at a very large energy scale (140 TeV), when the anomalous 
contribution starts to increase. 

Before drawing conclusions, it is necessary to comment on the other dependence, that 
on tan/3, which appears to be far less significant compared to Mi (or M' z ) and g%- This 
third parameter essentially has a (very limited) influence on the mass of the extra Z' and 
on the overall predictions. This is clearly shown in Figsj25] and [26], where we have varied 
both tan/3 and gs- Therefore, the mass of the extra Z' and the Stiickelberg mass may be 
taken to be essentially coincident, to a first approximation. 

• The bounds 

We conclude our analysis with two plots, depicted in Figsj27] and [28j which show the 
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Figure 26: Mass of the extra Z' gauge boson for different values of the Stiickelberg mass. 



variations of the bounds with the parameters qb and M\. In the first plots, shown in 
Figl27! we choose a large value of tan/? and we have varied both the Stiickelberg mass 
and the strength of the anomalous coupling. For a Stiickelberg mass around 1 TeV, the 
bound is around 5 TeV, for an anomalous coupling gs = 0.2. For a smaller value of 
gs = 0.1 the bound grows to 10 TeV. For a smaller value of gs = 10 -2 , the bound is 
around 100-120 TeV. This result is particularly interesting, because it should allow one 
to set limits on the Stiickelberg mass and the value of the anomalous couplings at the 
LHC in the near future. Smaller values of gn are tested in FigEH where the bound is 
shown to increase significantly as gn gets smaller. 



9 Conclusions 

We have analyzed the connection between the WZ term and the GS mechanism in the context 
first of simple models and then in a complete brane model, containing three extra anomalous 
U(l). We have shown that the WZ method of cancellation of the anomaly does not protect the 
theory from an excessive growth, which is bound to violate unitarity beyond a certain scale. We 
have also studied the connection between the two mechanisms, illustrating the corresponding 
differences. 

We have quantified the unitarity bound for several choices of the parameters of the theory. 
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Figure 27: Bound for different values of the Stiickelberg mass and tan j3 = 40. The bound grows as 
we reduce the anomalous coupling, and approaches the Standard Model behavior. For gs = 0.01 the 
model exhibits a bound around \fs max = 100 TeV (not shown). 

The significant dependences are those on the Stiickelberg mass and the coupling constants of 
the anomalous generators. We have also shown that the exchange of a physical axion lowers the 
cross section, but not significantly, whose growth remains essentially untamed at high energy 
We have shown that in these models one can identify a region that is SM-like, where some 
anomalous processes have a fast fall-off, from a second region, where the anomaly dominates. 
Clearly, more investigations are necessary in order to come out with more definitive predictions 
for the detection of anomalous interactions at the LHC, since our analysis has been confined 
to the parton level. The approach that we have suggested, the use of BIM amplitudes to 
search for unitarity violations at future colliders, can be a way in the search to differentiate 
between non-anomalous [30] and anomalous extra Z'. Our objective here has been to show 
that there is a systematic way to analyze the two mechanism for canceling the anomalies at 
the phenomenological level and that unitarity issues are important in order to characterize the 
region in which a certain theory starts to be dominated by the chiral anomaly. We hope to 
address these points in the near future in related work. 
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Figure 28: Bound for different values of tan/3, for Mi = 800 GeV. 
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10 Appendix. The cross section 

The total cross section is given by the sum of all the contributions shown in FigJTTI 



a MLSOM (s) 



a A (s) + cr B (s) + a c (s) + a D (s) + a E (s) + a F (s) 
+ ct A b(s) + a AE {s) + cr AF (s) + a BF (s) + cr EF (s). 



The interference term between Z exchange and x exchange, diagram (b) in Fi£ 
given by 



ctab(s 



1 

16tt 



l c i A ^ \ C 2 A SJ 



M 2 z (s - Ml) 



(120) 



TU), is 



121^ 



q J l / 

The following interference terms are vanishing 

M A M ] C + M C M\ = 0, M A M j D + M D M\ = 0, 

M B M 1 c + M c M B = 0, M b M ] d + M d M ] b = U. 



(122) 
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The interference term between the exchange of the Z boson and x exchange, diagram (e) in 
Fig. dUJ), gives 

vae(s) 



1024M|vr(s - M|) 



Y \ c " A ^ Y \°* A 



X 



L /' 



while the interference term between \ exchange and Z' boson exchange contributes with 



ctbe{s) 



16tt(s - M2) 2 

The other interference term in the cross section is given by 



^C (s,m/)c* 
. f 



^C (s,m q )c 



99 



L q 



<tcd(s 



> 9yy9n 



X-M 
99 



f 



16tt(s - M*) 2 *'"* 99 
so we obtain 

&cd(s) + (? BE {s) = 0. 
Other interference terms also vanish; in fact, we get 

m c m e + m e mIj = o, M d M e + M e M ] d = 0, 

M c M F + M F M 1 c = 0, M D M F + M F M ] D = 0. 
The interference term between Z exchange and Z' exchange takes the form 



(123) 

i 

(124) 

(125) 

(126) 

(127) 
(128) 



1 



oaf 



1024vr 



Y o c 2 A e,f 



MiMir (129) 



Y \ d " A ^ Y \ d ^f Y \ c i A ^' 

9 J L / J L q> J L /' 

The interference term between Z' exchange and x exchange of diagram (b) in Fig. (TTTT) is given 
by 



a B F(s 



16tt 



Y\ d " A ^ Y \ d 2 A ^f 



9na9 X 



M%,{s - Ml) 



(130) 



Finaly, the interference between Z' exchange and \ exchange, diagram (e) in Fig. (fTTj) . gives 



a EF [s\ 



1024M|,vr(s-M2) 



. r 



Y\ d " A ^ yV*^ 



(131) 
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